arXiv:1505.00567vl [hep-ph] 4 May 2015 


Preprint typeset in JHEP style - PAPER VERSION 


Next-to-Leading Order QCD Corrections to the Decay 
Width H Z 7 


Roberto Bonciani 

Dipartimento di Fisica, Universita di Roma “La Sapienza” and INFN Sezione di Roma, 

00185 Roma, Italy 

E-mail: roberto . bonciani@romal. infn. it 

Vittorio Del Duca 

Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland 
and INFN Laboratori Nazionali di Frascati, 00044 Frascati (Roma), Italy 
E-mail: delducavOitp . phys . ethz . ch.delducaSlnf . infn. it 

Hjalte Frellesvig 

Institute of Nuclear and Particle Physics, NCSR “Demokritos”, 

Agia Paraskevi, 15310, Greece 

E-mail: f rellesvigOinp. demokritos . gr 

Johannes M. Henn 

Institute for Advanced Study, Princeton, NJ 08540, USA 
E-mail: jmhenn@ias.edu 

Francesco Moriello 

Dipartimento di Fisica, Universita di Roma “La Sapienza” and INFN Sezione di Roma, 

00185 Roma, Italy 

E-mail: francesco .morielloSromal. infn. it 

Vladimir A. Smirnov 

Skobeltsyn Inst, of Nuclear Physics of Moscow State University, 119991 Moscow, Russia 
E-mail: smirnov@theory. sinp. msu. ru 

Abstract: We present the analytic calculation of the two-loop QCD corrections to the 
decay width of a Higgs boson into a photon and a Z boson. The calculation is carried 
out using integration-by-parts identities for the reduction to master integrals of the scalar 
integrals, in terms of which we express the amplitude. The calculation of the master 
integrals is performed using differential equations applied to a set of functions suitably 
chosen to be of uniform weight. The final result is expressed in terms of logarithms and 
polylogarithmic functions LD, Li 3 , LI 4 and LD^- 

Keywords: Higgs decay, Feynman diagrams, Multi-loop calculations. 




1. Introduction 


The discovery at the Large Hadron Collider (LHC) of the Higgs boson of the Standard 
Model (SM) [1, 2] calls for the investigation of the its properties, with the degree of agree¬ 
ment between the observed and the predicted behaviour being still an open question. The 
new boson decays into two photons or two electroweak W/Z bosons. It should also decay 
into a photon and a Z boson. If it is a SM Higgs boson with a mass of 125.1 GeV, the 
branching ratio is B(H —y Z 7 ) = 1.55 • ICC 3 , with an uncertainty of about 9% [3]. The 
H —y Z 7 decay is being searched by the CMS [4] and ATLAS [5] Collaborations at the 
LHC. 

The leading order evaluation of the SM H —y Z 7 decay was performed long ago [ 6 , 7]. 
As the Higgs boson has no electric charge, it does not couple directly to photons. Then 
the decay of a Higgs boson into a photon and a Z boson must be mediated at loop level 
by charged particles. At one loop, it is mediated by a heavy-quark loop [ 6 ] or a IF-boson 
loop [7], just as for the H —y 77 decay. H —y Z 7 may provide information on new physics, 
as a different H —> Z'y decay rate is expected if H is a non-SM scalar boson [ 8 , 9], or a 
composite state [10], or if different particles circulate in the loop [11, 12, 13]. 

The two-loop QCD corrections to the H —y Z'y decay were computed numerically 
in Ref. [14]. They correspond to two-loop corrections generated by a gluonic exchange 
within the heavy-quark loop. In this paper, we present an analytic calculation of the QCD 
corrections to H —y Z'y. 

The paper is organised as follows: in Sect. 2 , we recall the definition and the explicit 
expressions of the amplitude for the decay H —y Z'y at leading order; in Sect. 3, we describe 
how we perform the analytic computation of the NLO QCD corrections to the decay width; 
in Sect. 4, we provide numerical results as a function of the Higgs mass; in Sect. 5, we draw 
our conclusions. The appendices contain details on the projection operators used to define 
the decay amplitude (App. A), the master integrals (App. B), the matrices of the system 
of differential equations for the master integrals (App. C) and the analytic properties of 
the functions occurring in the master integrals (App. D). 

2. H — y Z'y at Leading Order 

2.1 The amplitude for the decay H —y Z 7 

At leading order the SM Higgs boson decays into a photon and a Z boson via either a 
heavy-fermion loop or a IF-boson loop. The corresponding Feynman diagrams are shown 
in Fig. 1. Let us label the Z-boson momentum as p\ and the photon momentum as p 2 . 
The general Lorentz structure T^ u of the amplitude, 

M = T^£^ Pl )£ u {p 2 ), (2.1) 

for the decay of a Higgs boson into a Z boson and a photon with polarization vectors S/J/pi) 
and £ v {p 2 ) respectively, is 

T ,IV = p'M Ti + pW 2 T 2 + T 3 + T 4 + 5^ T 5 + P p 2 ,T 6 , (2.2) 
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(a) 


(b) 


(c) 


Figure 1: Feynman diagrams for the decay process H — > Z 7 at leading order. 


where p 2 = Pi ~Po an d Pi = P 2 = (Pi + P 2) 2 = — m 2 H . mp is the Higgs mass 

and rriz is the Z-boson mass. The coefficients T) are functions of the Mandelstam invariants 
of the problem under consideration and can be extracted using projector operators Pj 11 ' 
such that T{ = Pi tlu T llu . The projector operators are collected in Appendix A. 

Requiring the photon gauge invariance, T^ v p 2v = 0, in Eq. (2.2), we obtain 

T\ = 0, 

2 2 

m _ m T r — 171 % 

T5 = — pi ■ P2T1 = ---— T4 • (2.3) 

Furthermore, T 2 , T 3 and Tq* do not contribute to the squared amplitude |Ad| 2 . Thus, up 
to contributions which vanish in |Ad| 2 , the form factor can be written as, 

T^ = (p^-pi-p 2 6n^. (2.4) 


2.2 Leading Order Contribution 


The width for the decay of a Higgs boson into a photon and a Z boson can be cast in the 
following form, 


r^z 7 


G F a 2 (■ m 2 H - m |) 3 2 

64\/27r 3 miy 


(2.5) 


where a is the fine structure constant, Gp is the Fermi constant and where we introduced 
the function T related to the form factor X 4 by the following equation, 


T = 


167T 2 myy 

„ 3 S 2 

y s w 


t 4 . 


( 2 . 6 ) 


In Eq. (2.6), g is the weak coupling constant, s\y = sin 9w is the sine of the weak mixing 
angle, and mw is the IT-boson mass. 

‘Note that Te receives contribution only from the axial-vector part of the ZQQ vertex. However, since 
the Higgs particle is a C-even state, and the photon is C-odd, only the C-odd coupling of the Z (i.e. the 
vector coupling and not the axial-vector one) contributes to the decay width. This implies that T§ does not 
contribute to the general form factor (2.2). In fact, when we add together the diagrams with an opposite 
flow of the fermionic arrow, the contribution to Te changes in sign, in such a way that in the sum Te = 0. 


- 2 - 








The function T can be expanded in powers of the coupling constants, starting with 
the one-loop contribution. At leading order we have, 

f rp3 


= £k j-w +NcQ 


Qq s w 


J a ? 


(2.7) 


s\y swcw 

where N c is the number of colors, q = t,b labels the type of heavy quark, either a top or 
a bottom, circulating in the loop, Q q is the heavy-quark electric charge in units of e, T 2 
is the third component of the heavy-quark isospin, and cw = cos 9w is the cosine of the 
weak mixing angle. The actual expressions for and J- q 11 ^ are, 

(m 2 H m 2 z + 2m 2 z my V - 2m 2 H rrtf v - 12 m^) 


tt (^0 _ 


+ 


+ 


m w( m2 H ) 

/OOO /tO /IO 0/1\ 

(2 mjjm z myy — 2m|m^ — m z m H — 12m z rrtf v ) 


m 2 H (m 2 H — m 2 z ) 2 
(4 m z — 2m 2 z m 2 H + 12 — 12 m^m^ — 24 m^y) 


I{xw,Dw) 


( m 2 H — m 2 z ) 2 


J(xw,Uw), 


( 2 . 8 ) 


8 m 2 


— 

q (' m 2 H — m 2 z ) ' ( m 2 H — m 2 z ) 


+ 


8 m 2 m 2 7 

-2T2 1 ( x g,y q ) + 


4m 2 (4m 2 — m 2 H + m 2 z ) 


( m ‘H ~ m z) 2 


with 


4 m 2 


1 (*/. Vf ) = \ 1 - —2" lo 6(®/) - \ 1 - 


Am 2 , 


m 


H 

J(xf,yf) = ^ log 2 (x f ) - ^ log 2 (y/) 
where the variables x/ and yf are defined through, 

,2 (! ~ x ff 


rn 


!og (yf) 


H 


m 2 H = —m 


r 


2 2 (-*■ VfY 

m z = — m/- £ — 


J(x q ,y q ),(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


*/ yf 

with / = W,q and q = t,b. 

Let us consider real values of m t , and m#. In the region 0 < rrifj < 2m /, Xf 

has an imaginary part, with unit modulus and phase between 0 and ir, 


Xf = exp < i arctan 


yJm 2 H (Am 2 — m 2 

2m 2 , — m 2 rr 


(2.13) 


l' j I / L 

In the region mu > 2m/, we have — 1 < Xf < 0; we define Xf = —x'f + i0 with 


x f = 


m 2 H — Jm 2 H — Am 2 


\Jm 2 H + \Jm 2 H — Am 2 
Analogously, in the case 0 < mz < 2mf the variables yf are on the unit circle 


(2.14) 


yf = exp < i arctan 


m|(4m 2 — m 2 z ) 
2 m 2 — m 2 z 


1 


(2.15) 


- 3 - 




(d) (e) (f) 

Figure 2: Feynman diagrams for the NLO QCD corrections to the decay process H —> Z 7 . 
Diagrams with the reversed direction of the fermionic arrow are not shown. We can easily consider 
the N% contribution (one top-quark loop coupled to the Higgs and the other correcting the Z or 7 
propagator). 


while for mz > 2 m f we define yj = —y'j + iO with 



3. NLO QCD Corrections 


(2.16) 


The diagrams involved in the calculation of the NLO QCD corrections to the decay width 
of a Higgs boson into a photon and a Z boson are shown in Fig. 2. Their contribution to 
the form factors can be extracted using the projectors defined in Appendix A. Expanding 
in the strong coupling constant we have, 

T = F {11) + —J^ 2Z) + ..., (3.1) 


(2Z) 

At this order in as, the bare form factor, Jq , is UV divergent and needs to be renormal¬ 
ized. The only renormalization required is the heavy-quark mass renormalization, in the 
fermionic propagators and in the coupling of the Higgs boson to the heavy-quark pair. We 
choose to perform the mass renormalization in the on-shell (OS) scheme [15], 


0m OS l e i m qi 2 


a S n( , 
— m q —G(e) 

7 T 



C F (3 — 2e) 
T e (1 — 2 e) ’ 


(3.2) 
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Figure 3: Counterterm diagrams involved in the heavy-quark mass renormalization. 


where /i is the renormalization scale, Cp = (TV 2 — l)/(2iV c ) with N c the number of colors, 
(7(e) = (47r) e r(l + e) and e = (4 — d)/ 2, with d the space-time dimension. Indicating with 
the renormalized form factor, we have 

F {2l) = F^ l) + CT . (3.3) 


The contributions to the counterterm come from the diagrams shown in Fig. 3. 
Retaining only terms of 0(as), Eq. (2.5) can be written as, 


FH^Z'y = 


Gfj,oi 2 


64\/27r 3 miy 


K(JF( i; ))gfi(J-( 2 0) + cj(jr( 1 0)cj(jr( a )) J 


H 
etc 

+2 — 

7r 


(3.4) 


— + &qcd) , 


where we defined 

n a s ^{F^)^{F^) + ^(F^)Q(F^) 
QCD ~ 7 T 3?(^)) 2 + ^( jr ( tO ) 2 

as of the NLO QCD corrections with respect to the LO contribution. 


(3.5) 


3.1 Calculation of the master integrals 

jt(2Z) 

is expressed in terms of a large number of scalar integrals that are individually ultra¬ 
violet divergent. To deal with these divergences, we perform the integrals in dimensional 
regularization [16, 17, 18, 19, 20, 21]. The dimensionally regularized scalar integrals are not 
all independent. The reduction to a set of independent integrals, called master integrals, is 
carried out using two different computer programs, FIRE [22, 23, 24] and Reduze [25, 26], 
that implement in an automatic way the solution of the linear system of integration-by¬ 
parts identities [27, 28] which relate the Feynman integrals. 
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p2 Q p 

(Mx) (M 2 -M 3 ) (M4-M5) (M 6 -M 7 ) 











Figure 4: Set of master integrals for the NLO QCD corrections to the decay H —»• Z 7 . The master 
integrals M2, M4, Mg, Mg are functions of p 2 = —m 2 H , while the master integrals M3, M5, M7, 
M 10 are functions of p 2 = —m 2 z . M 15 and Mig have the numerator explicitly written. A dot on a 
propagator indicates that the propagator is raised to power 2. Two dots means that the propagator 
is raised to power 3. See Appendix B. 
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We calculate the master integrals using the Differential Equations Method [29, 30, 31, 
32, 33, 34, 35], following recent developments [36, 37, 38, 39, 40] (see also [41, 42, 43, 44, 
45, 46] and [47, 48] for further studies of the method). 

The system of linear differential equations is cast, via a suitable basis choice for the 
master integrals, in the canonical form [36], 

df = e dA f . (3.6) 


The e dependence is completely factorized from the matrix dA. The latter depends only 
on the dimensionless variables x q and y q , defined in Eq. (2.12). 

f is a vector of 28 functions, fi(x q ,y q ,e), defined in terms of the integrals drawn in 
Fig. 4, 


h 

h 

h 

h 

fh 

k 

h 

h 

h 

ho 


fu 

fl2 

/l3 


fu 

/l5 

/l6 

fl7 


16 e 2 Mi , 

16 e 2 rrijf Mu , 

8 e 2 sJ-m 2 H {4m 2 - m 2 H ) (M 4 + 2M 6 ) , 
16 e 2 m| M 5 , 

8 e 2 \J—m 2 z {4m 2 - m 2 z ) (M 5 + 2 M 7 ) , 
16 e 2 yj—m 2 z (4m. 2 — m|)Af 3 , 

16 e 2 \J—m 2 H (4m 2 — m^)M 2 , 

16 e 3 (m# - m|) M, 6 , 


16 e 2 m 2 (m# — m 2 z ) Mi 7 , 



1 - 


4m^ 




[m^(4m 2 Mi 7 - M 4 ) + 4m|(m 2 M 17 - M 18 ) 


—2 e (m^ + m'l) M^] , 

16 e 3 (m^ - m|) M 13 , 

16 e 2 m 2 (m^ — ru.|) M 44 , 

1 r 4m 2 


4 e 


2rri 2 — m 


1 - 


m 




2 4m 2 m 2 H - m 2 H m 2 z + m% 
m z --- 5 - 5 - m 5 


—4 m 2 H 


m 2 m 2 H — m 2 H m 2 z + m z 

2 2 

m H — m z 


m H ~ m 2 


M 15 - 4 m 2 (m 2 H - m|)m| M 14 


+2e(2 m 2 m 2 H + m 2 H m z - m|)M 13 


16 e 3 (m^ - m|) M 12 , 

16 e 2 m| (4m 2 - m|) M 10 , 

16 e 2 yJ—m 2 H (4m 2 — m 2 H ) \J—m 2 z (4m 2 — m|) Mu 
16 e 3 (rn 2 H - m|) \J-m 2 z (4m 2 q - m 2 z ) M 20 , 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


(3.16) 

(3.17) 

(3.18) 


(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 


- 7 - 



/is = 16 e 4 (rnjf - m z ) M 2 1, (3.24) 

fig = 16 e 3 (rn 2 H - m z ) -m 2 z {Am 2 - m|) Af 22 , (3.25) 

/ 20 = 16 e 4 (m 2 H - m|) M 25 , (3.26) 

/ 2 1 = -16e 3 m| (rr 2 H - m 2 z ) -m 2 z {4m 2 - m|) Af 26 , (3.27) 

/ 22 = 16 e 3 (m 2 H - m|) -m 2 H {4m 2 - m 2 H ) M 27 , (3.28) 

/ 23 = 16 e 2 [2{m 2 H m 2 z - 2m 2 m 2 H - 2m 2 m 2 z )Mu + m 2 (m 2 H - m 2 z ) 2 M 28 

+ e(m 2 H - m z )(m z M 26 - M 27 )\ , (3.29) 

/ 24 = 16 e 4 (m 2 H - mn 2 z ) M 23 , (3.30) 

/25 = 16 e 3 (m^ - m|) -m 2 H (4m 2 - m 2 H ) M 24 , (3.31) 

/ 26 = 16 e 3 (rn 2 H - rn z ) M 8 , (3.32) 

/ 27 = 16 e 2 m 2 H (4m 2 q - m 2 H ) Mg , (3.33) 

/ 28 = 16 e 3 (m^ - m|) iJ-m 2 H (4m 2 - m 2 H ) M ig . (3.34) 


The functions fi(x q ,y qi e) are chosen in such a way to be pure and of uniform weight, in 
the sense of [36]. 

The explicit definition of the integrals M \, • • • , M 28 is given in Appendix B ' . To set 
the normalization, we define M\ as the following integral, 


Mi 


= / V D ki V D k 2 


(kf + m 2 ) 2 ( k 2 + m 2 ) 2 


where the integration measure is such that, 


1 

16?’ 


(3.35) 


d D ki = 4tt 2 -T(1 + e)(^ ) V D ki . 

\m q J 


With this normalization, 


/i = l- 


The matrix dA is a differential depending on x q and y q , 


(3.36) 


(3.37) 


dA = Si dlogXq + S 2 d log (1 — x q ) + §3 d log (1 + x q ) + § 4 d log y q + §5 d log (1 - y q ) 
+§6 d log (1 + y q )+S 7 d log (x q - y q ) +§ 8 d log (1 - x q y q ) 

+§9 d log (1 - Xq- x q y q + x 2 ) + §i 0 d log (1 -y q - X q y q + y 2 ) 

+§11 dlog (x q -y q + x q y q - x 2 y q ) + §i 2 dlog (x q -y q - x q y q + x q y 2 ). (3.38) 


The sparse matrices §,; are purely numerical and they are collected in Appendix C. 

To fully describe the f, we need to complement the differential equations with boundary 
conditions. We can, in principle, choose any kinematic point (x q ,y q ). However, it is well 
known that often, (spurious) singularities of the differential equations allow one to fix the 


'The calculation of some of the master integrals in Fig. 4 was already performed in Ref. [49]. 



boundary condition without calculation, using the physical insight that certain limits must 
be non-singular. This is also the case here. We notice that the integrals defining the primary 
basis (Fig. 4) are regular in the point s = m 2 z = 0. Moreover the rational prefactors of the 
combinations given in Eqs. (3.7-3.34) vanish in the same limit. Therefore, all the integrals 
fi(x q , y q ), i = 2,..., 28, vanish in this limit. The only exception is represented by /i, which 
is identically equal to 1. We can write the full set of boundary conditions in the compact 
form, 

/i(l, 1) = <y M . (3.39) 

The system (3.6), together with the boundary condition (3.39), makes it obvious that the 
solution, i.e. the functions fi(x q , y q ), have a number of desirable properties. At any order in 
the expansion in e, they are given by iterated integrals [50] over the one-form dA. Defining 
the weight of an iterated integral as the number of integrations, we see that at order e fc , 
f is given by a Q-linear combination of iterated integrals of weight k. Such functions are 
referred to as pure functions of uniform weight. 

The basis choice for f leading to this form was achieved using the ideas outlined in 
Ref. [36]. Specifically, generalized unitarity cuts in four dimensions were used to project 
onto subsets of the differential equations. This typically leads to an answer close to the 
canonical form, where unwanted terms e.g. due to integrals vanishing on the cuts can be 
easily removed algorithmically, see e.g. Refs. [41, 42], A canonical basis for very similar 
integrals was found in Ref. [41], and using the results from that paper the present basis 
choice was rather straightforward. 

The one-form dA characterizes the type of iterated integrals that are needed. In 
particular, given the definition dA = Yl\=\ §idlog(aj), the individual a* are called the 
letters, and iterated integrals correspond to words in those letters. The set {«*} is called 
alphabet. 

For the specific one-form dA, the system in Eq. (3.6) can be solved in terms of a 
subset of generalized (or Goncharov) polylogarithms (GPL) [51, 52, 53], defined as iterated 
integrations over a set of basic polynomials {(x — ai),..., (x — a n )}, 


G(0 n ,x) = — log n x, G(ai,a 2 , ...,a n ,x) = / G(a 2 ,a n , t ), (3.40) 

n\ J 0 t - ai 

for a certain set of a*. 

We choose to represent the GPLs of two variables, x q and y q , as functions of argument 
x q and weights depending on y q . The weights for the GPLs functions of x q are determined 
by the following set, 

jx q , 1- x q ,l + x q ,x y - y q , x y - 1 /y q ,x q - R L | , (3.41) 

where the roots Rj are defined by the following expressions, 


Roi 


R02 


Vq 

vl 


Vq 

- y q +1 ’ 

~ Vq + 1 
Vq 


(3.42) 
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(3.43) 



Rii ~\ + 2j q ~ ^y i+2yq - 3y ^ (3 - 44) 

Rl2= 2 + 2V q + 2^J l + 2yq ~ 3y2q ’ (3 - 45) 

-^21 = - + ~y q — - \J —3 + 2y q + , (3.46) 

^22 = - + + ~\J—3 + 2y q + y% . (3.47) 


We also have GPLs of the variable y q . The weights of these GPLs are determined by 
the following set, 


where 



y q , 1 + y q , y q - G y q - C, y q - i, y q + ij , 

(3.48) 

1 — iy/ 3 _ 1 + i\J 3 

2 ’ C ~ 2 ’ 

(3.49) 


are sixth roots of the unity. 

The solution of the differential equations can be found expanding each of the canonical 
master integrals in e, 


fi 


if + PWifW 


(3.50) 


We can identify the terms on each side of Eq. (3.6) which multiply the same power of e, 
i.e. 

df\ k) = dAijff - 1} . (3.51) 

This allows us to solve Eq. (3.6) recursively order by order. As the values of /are 

constants given entirely by the boundary conditions, we shall start by finding f- V> (x q , y q ). 

This may be done by integrating the right-hand side of Eq. (3.51) from the boundary point 

(1,1) to a general point (x q ,y q ), and we choose to use a stepwise path passing through 

(x q , 1). Using the boundary conditions to fix the integration constant, gives the values for 

f,- 1 ) ■ Once this is done for each of the canonical master integrals, the procedure may be 
( 2 ) 

repeated for f- ' and the other orders, yielding the full result for /) up to any desired order. 

We checked all the expressions of the master integrals numerically against the computer 
program Fiesta [54, 55, 56] finding complete agreement. 

The analytic expressions of the functions /) are collected in an ancillary file of the 
ai'Xiv submission of the present paper. 

A comment is due regarding analytic continuation and analytic boundary values in 
different regions. In the following subsections, we explain how we obtained a formula 
valid in all kinematic regions. For completeness, here we wish to comment that another 
possibility is to solve the differential equations directly starting from a boundary value 
in the region of interest. The boundary value we gave in Eq. (3.39) at x q = y q = 1 is 
in a non-physical region, and other regions can be reached by analytic continuation. In 
particular, when analytically continuing to negative values of x q or y q , the d\ogx q and 
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d\ogy q terms in dA are responsible for imaginary parts. As an example, let us analyze the 
boundary point (x q , y q ) = (—1, —1). Unlike {x q ,y q ) = (1,1), this is a singular point, so that 
it has to be approached with care. Since x q = y q is nonsingular, we can set x q = y q = z, 
for z > {) in which case the alphabet becomes {dlog z, dlog(l — z),c21og(l + z)}. When 
analytically continuing to negative values of z, we have to take care of the imaginary parts 
due to dlogz, keeping in mind that z has a small imaginary part ?'0. In this way, one can 
give an analytic boundary value as z —>• — 1 , in terms of powers of log(l + z) and a set of 
constants. It follows from the reduced alphabet {dlog z, c?log(l — z),d log(l + z)} that the 
only constants required are Euler sums. Up to weight four, one has iir, log 2, £ 3 , Li 4 (l/ 2 ), 
and products thereof. 

3.2 An alternative choice of the functional basis 

The analytic expressions of the master integrals in terms of GPLs provide all the necessary 
features of an analytic formula: possibility to expand the formulae in particular regions 
of the phase space, flexibility with respect to the input physical parameters, available 
routines for their numerical evaluation [57]. However, for the sake of a faster and more 
stable numerical evaluation, we describe in this section how to represent the result in terms 
of a different functional basis. 

Up to weight four, Goncharov polylogarithms can be rewritten in terms of the following 
functions (see [58] and references therein), 

log(aq) = G(0,aq), Li n (aq) = -G(0„_i, 1, aq), Li 22 (aq, x 2 ) = G ( 0, —, 0, , 1 ) , 

V X 2 X 1 X 2 J 

(3.52) 

with n = (2,3,4), and where x\ and x 2 are rational functions of x q and y q . Li n (xi) has a 
branch cut for x\ > 1 , while Li 2 ) 2 (aq,x 2 ) has a branch cut whenever x\ > 1 or xix 2 > 1 . 
With this choice, the functions given in Eq. (3.52) can be directly evaluated using the 
numerical routines of [57]. 

In order to find an expression in terms of these functions, the concept of the symbol 
[51, 59, 58] of an iterated integral is very useful. The symbol corresponds to the integration 
kernels defining the iterated integrals. It is completely manifest in our differential equations 
approach. It is possible to use symbol-based ideas to rewrite the master integrals in terms 
of a minimal function basis, see e.g. [60]. Moreover, we can recover the information about 
the terms in the kernel of the symbol using the coproduct map [61]. 

However, it is also possible to proceed in a more direct way, using the knowledge of the 
differential equations, proceeding in the following algorithmic steps. First, one generates a 
list of function arguments as monomials in the letters appearing in our function alphabet 
(i.e. the arguments of the logarithms appearing in the one-form dA.) For the classical 
polylogarithms Li n (xi), one requires that 1 — x\ factorizes over the letters appearing in 
the alphabet. (A caveat is that in principle, “spurious” letters might be needed [60], but 
this was not the case here.) For Li 2 ) 2 (aq,x 2 ), the condition is that 1 — x±, 1 — x 2 ,1 — xix 2 
factorize over the alphabet. Second, for each weight k, one chooses a maximal set of linearly 
independent functions for the alphabet (the linear independence can be verified using the 
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symbol). By construction, the differential equation at oder e k can then be solved in terms 
of this set of functions. 

The Li expressions for the master integrals are provided in an ancillary file and are 
defined for the relevant physical region s > rn 2 z , which in terms of the new variables implies 
y < x or arg(y) < arg(x), where 0 < arg(x) < ir is understood. 

When solving the differential equations in terms of the above functions, there is a 
certain freedom in the choice of the set of function arguments. The choice is usually 
guided by trying to make certain properties of the answer manifest, such as simple branch 
cut properties. For instance, one could require that the functions are real valued in the 
physical region, so that the imaginary parts of the master integrals become explicit. We 
found that this was possible in the regions of real x q ,y q . However, for x q or y q complex, 
individual functions develop imaginary parts. 

In order to obtain an expression valid over the entire physical domain, a linear com¬ 
bination with rational coefficients of the Li functions is not sufficient. The reason is that 
some of the functions in Eq. (3.52) have discontinuities in the physical region, and in order 
to ensure the analyticity of the master integrals it is necessary to introduce discontinuous 
functions (Heaviside theta functions) that cancel branch cut discontinuities of the new ba¬ 
sis. The same issue was discussed in Ref. [60]. In Appendix D we briefly introduce the 
issue, and we provide a method to find an expression for the master integrals in terms of 
Li functions valid over the entire physical region. 

4. Numerical Results 

In order to assess the impact of the NLO QCD corrections, it is convenient to recall 
the relative contributions of the top, bottom and electroweak loops to the leading order 
width. For mn = 125.1 GeV, mt = 173.34 GeV, = 4.6 GeV, m\y = 80.398 GeV, 
rn z = 91.1876 GeV, = 0.23149, a = 1/128, G F = 1.16637 • 1(T 5 , the QCD width is 
0.02 KeV at leading order. In Table 1, we report the relative contributions of the top- 
quark loop, the bottom-quark loop and the top-bottom interference to the leading order 
QCD width. We note that the bottom loop contributes only one per rnille to the QCD 
width, which is almost entirely given by the top loop, with a sizeable destructive top- 
bottom interference. On the other hand, the full width is 6.67 KeV at leading order. In 
Table 2, we report the relative contributions of the QCD (i.e. top and bottom) loops, 
the IT-boson loop, and the QCD-lT-boson interference to the leading-order width. We 
see that the QCD loops contribute only three per rnille to the width at leading order; the 
IF-boson loop accounts for the bulk of the width, up to a large destructive QCD-lT-boson 
interference. 

The NLO QCD corrections to the decay width of a Higgs boson into a Z boson and a 
photon were computed numerically in Ref. [14] and turn out to be quite mild. In Fig. 5, 
we plot the decay width as a function of the Higgs mass. T H^-z-y includes the NLO 
contributions of the diagrams with a top-quark loop and a bottom-quark loop. In Fig. 6, 
we plot the two-loop QCD corrections to the decay width. In particular, we plot the 
contribution of the top loop alone, the bottom loop alone, and the complete correction. 
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r { QCD (KeV) 

p(i0 /p(U) 

1 t / 1 QCD 

■p(10 /T"'(10 

L b / l qcd 

Interf ™/ r QCD 

0.02 

1.052 

1 •10~ 3 

- 0.053 


Table 1: Values of the QCD leading-order width, and of the relative contributions of the top loop, 
the bottom loop and of the interference between top and bottom loops, at run = 125.1 GeV. 


T^ (KeV) 

r QCr>/ r(U) 

7^ 

t-H 

Interf QGZ?-uV r{li) 

6.67 

3 • 10" 3 

1.112 

- 0.115 


Table 2: Values of the leading-order width, and of the relative contributions of the QCD loops, 
the TV-boson loop and of the interference between the QCD and the IV-boson loops, at uih = 
125.1 GeV. 



rriH (GeV) 


Figure 5: Tu^z-y including LO and NLO QCD contributions. 


For small values of the Higgs mass, the corrections due to the top loop and those due to the 
bottom have opposite sign, with the top correction being dominant. In Table 3, we quote 
the figures for mu = 125.1 GeV and as(m 2 H ) = 0.115: the two-loop QCD corrections 
amount to 0.22% of the leading order width, with the top loop adding a 0.3% and the 
bottom loop subtracting a 0.08%. The two-loop top corrections is in agreement with the 
corresponding plot of Ref. [14] and with the results in Ref. [62] 

As we recalled above, at one loop the top and bottom loops contribute 0.3% of the 
leading order width. Since the two-loop QCD contribution is about 75% of the one-loop 
QCD contribution, that implies a bad convergence of the QCD part of the width. In 
principle, higher QCD orders would be necessary to stabilise the QCD perturbative series. 
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rriH (GeV) 

Figure 6 : The size of the NLO QCD corrections with respect to the leading order contribution, 
5qcd i as defined in Eq. (3.5). In this figure we plot the contribution of the top loop only, the 
bottom loop only, and the top and bottom loops. 


T (KeV) 

Sqcd (top) 

Sqcd (bottom) 

Sqcd 

6.68 

3 • 10" 3 

-0.8 • 10~ 3 

2.2 • 10" 3 


Table 3: Values of the NLO QCD width, and of the NLO QCD corrections Sqcd with respect to 
the leading order contribution, at mn = 125.1 GeV. 

However, up to NLO accuracy, the QCD contribution is numerically tiny. Thus, it would 
be more relevant to compute the two-loop electroweak corrections, which we expect to be 
larger than the QCD ones. 

5. Conclusions 

In this paper we present the analytic calculation of the NLO QCD corrections to the width 
of a Higgs boson decaying to a Z boson and a photon. These corrections were computed 
numerically in [14], where the authors studied the size of the corrections that come from 
the two-loop Feynman diagrams with a top-quark loop, with respect to the leading order 
contribution. In our paper we include also the contribution due to the two-loop diagrams 
with a bottom-quark loop. For small values of the Higgs mass, the corrections due to the 
top loop and those due to the bottom have opposite sign, with the top correction being 
dominant. For my = 125.1 GeV they amount to 0.22% of the leading order width. 

The calculation was carried out using integration-by-parts identities for the reduction 
to master integrals of the dimensionally regularized scalar integrals, in terms of which we 
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expressed the amplitude. The calculation of the master integrals was performed using the 
Differential Equations Method applied to a set of functions suitably chosen to be of uniform 
weight. 

The solution is expressed in terms of logarithms and polylogarithmic functions Li n , 
with n = 2, 3,4, and Li 2 , 2 - Their arguments are rational expressions in the two dimension¬ 
less variables x q , y q , written in terms of the three mass scales of the problem: rrifj, mz 
and m q . The numerical evaluation of the polylogarithmic functions is done using existing 
numerical C++ routines. 

Finally, we recall that for rriH = 125.1 GeV, at leading order the top and the bottom 
loops account for about 0.3% of the width, the overwhelming contribution being yielded 
by the electroweak loop. Having found that the two-loop QCD corrections, which correct 
the heavy-quark loop, equal about 0.22% of the leading order width, it is reasonable to 
expect that the two-loop electroweak corrections, which correct the electroweak loop, will 
yield a larger contribution +. 
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*The Z boson may decay leptonically - the final state of H — > Z 7 is a lepton pair, l~l + , and a photon - 
or hadronically - the final state is formed by two jets and a photon. However, the two-jet plus photon final 
state is overwhelmed by the QCD background. Conversely, because of the smaller pp —> l~l + 7 background, 
the decay H —> Zy —> l~l + 7 is much cleaner. In order to compare to the pp —> l~l + 7 background, it 
is different to consider the resonant production H Z 7 — > l~l + 7 or the off-resonant one H —¥ l~l + 7 , 
for which no on-shcll Z-boson is produced in the intermediate state. Estimates of the difference between 
the resonant and the off-resonant widths vary [63, 64, 65, 66 ], however with realistic experimental cuts, 
the difference between the resonant and off-resonant leading-order widths for electron or muon pairs is of 
the order of a few per cent [67]. As regards the two-loop QCD corrections we computed, the off-resonant 
diagrams are not present. Thus, with little effort - just including the Z-boson decay into two leptons - 
we could have obtained the QCD corrections to the decay H —> l~l + 7 . However, as regards the two-loop 
electroweak corrections, there are also off-resonant diagrams. Therefore, expecting that the still unknown 
two-loop electroweak corrections may be larger than the two-loop QCD corrections, we think it is best to 
postpone an analysis of the H —¥ l~l + 7 decay to the computation of the two-loop electroweak corrections. 
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A. Projection operators 

The projector operators to be used in Eq. (2.2) are 


{my - m 2 z ) 


LL 1/ 

2Y2 P2P2 


16 {d — 1 )m z 


(m 2 H - m\Y PlPl + {d - 2){m? H - m|) 4 P * P2 

8 ( d ~ W , _#*_*!_ 4m l *«/ 

{d — 2){m 2 H — m 2 z ) 3 12 21 (d — 2){m? H — m 2 z ) 2 

8{d - l)m 2 z ftv _4_ // v 

{d— 2){m 2 H — m 2 z ) 3 d>2 ^ >2 {d — 2 ){m 2 H — m|) 2 ^ 1 ^ 2 

_L_ ~ I_^_ fiH” 

{d— 2){m 2 H — m 2 z ) 2 2 1 (d — 2){m 2 H — m 2 z ) 

8(d-l)m? z „ 4{d-l) 

{d — 2){m 2 H — m 2 z ) 3 2 2 (ci — 2){m 2 H — m 2 z ) 2 1 2 

(d — 2){m 2 H — m 2 z ) 2 P2Pl (ci — 2){m 2 H — m 2 z ) 

Am ^ 2 

-7^7-^-Z 2 W P 2 P 2 + 7T K^2 + 


(d — 2){m 2 H — m z ) 2 

4_i_ 8^ 

(d~ 2) ’ 

4 


(d~2)(m 2 H -m|) 


{d —2 ){d — 8){m 2 H — m 2 z ) 2 


e^P a PlpP2( 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 


B. Primary master integrals 

In this appendix we give the explicit expression of the “primary” master integrals shown 
in Fig. 4. 

We define the following set of seven denominators, 

Di = kf + m 2 q , 

D 2 = k 2 , 

D 3 = (fa + k 2 ) 2 + m 2 , 

D 4 = {pi -kif + m 2 , 

D 5 = (p 2 + h ) 2 + m 2 , 

Dq = (p\ - k x - k 2 ) 2 + m 2 , 

D 7 = {p 2 + ki + k 2 ) 2 + m 2 q . 

The integrals of Fig. 4 are then, 

M l =fV>k l V>k,^ SI = (B.8) 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 
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M 3 
M 5 
M 7 
Mg 
Mu 
M13 
Mis 
Mi 7 
Mi g 
M 2 i 
M 2 3 

m 25 

M 2 7 


1 

D 2 2 D 2 D 4 


V D k x V°k 2 


S v ° k ^w^' 

S vDk ^ Dk "-w^' 


V D h V D k 2 


d 2 d\ d 4 d 5 ’ 

/u’l + m 2 


f n n k i + m o 

J° Dk ^n^o(o e ' 

[ V D h V D k 2 1 
J D\ D 4 D 3 Dq D 

[ V D ki V D k 2 --- 

J 2 D 4 D 2 D 3 D 5 D 

[ V D kx V D k 2 --- 

J 2 D 2 D 3 D 4 D 5 D 

[ V D k i V°k 2 --- 

J D 2 D 3 D 4 D 5 D 

f V D h V D k 2 --- ir— 

J 2 d 2 d 3 d 4 d 2 d 


m 4 

m 6 

Ms 
Mio 
Mi 2 
Mu 
Mi 6 

Mis 

M 20 

m 22 

M 2 4 

m 26 

m 28 


= J V D h V D k 2 
= J V D ki V°k 2 
= J V D ki V D k 2 
= J V D k\ V D k 2 

= j V D k\ V D k 2 
= j V D ki V°k 2 

= J V D kx V D k 2 

= j V D kx V D k: 2 

= J V D k\ V D k 2 
= J V D k\ V°k 2 
= J V D ki V D k 2 


1 

d-2 b\ d 2 ’ 

1 

Dl D\ d 6 ’ 

1 

/) D 2 /) T Z?5 
1 

Z>f d| d 4 Dq ’ 

i 

A A), D 5 D| ’ 

1 

D 2 I)'\ d 4 d 3 

1 

D 2 D3 D 2 Dg ’ 

(g2 + fcl + fc2) 2 

d 2 d\d 2 d % ’ 

1 

Dl D 2 D 4 D 5 D 6 ’ 
1 

Dl d 2 D 3 D 5 Df ’ 
1 

D 2 D3DID5D7 ’ 


= / V D ki V D k 2 


= / V D k\ V D k 2 


D 2 Do D4 D5 Dg 


D 2 Do D 4 D|D 6 ' 


(B.9) 

(B.10) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 

(B.19) 

(B.20) 

(B.21) 


C. Matrices for the System of Differential Equations 


In this appendix we collect the matrices §i-§i 2 defined in Eq. (3.38). 
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D. Solution of the differential equations in terms of Li functions 

By definition the solution of the differential equations can be valid only in a region were the 
differential equations are non-singular. When a singular point of the differential equations 
is reached the solution might develop a branch cut, depending on the functional basis used 
to represent the solution. However in a physical process the master integrals are analytic 
in the physical region and those branch cuts must be spurious. In order to ensure the right 
analytic properties of the master integrals, we can replace the discontinuous functions with 
different, analytic branches. 

Solving the differential equations in terms of Li functions, we can in principle find a 
basis without branch cuts in the physical region, so that the solution has the right analytic 
structure and no spurious branch cuts are present. However in our case the set of analytic 
functions in the physical region is too small to represent the master integrals, and we are 
forced to employ also discontinuous functions. For simplicity they are chosen to be the 
following two logarithms and five Li3, 

log(xq - y q ), log (-y q x q + y q x q -y q + x q ) , Li 3 (a;) i = 1,..., 5 . (D.l) 

where the cq’s are 

x q y q - x q y q + x q -y q - x q y 2 + x q y q -x q + y q -x q y\ + x q y q -x q + y q 
(x q + l)(l-y q ) 2 ’ (1 - x q )(x q + l)y q ’ (x q + l)(y q -x q )(l - x q y q ) ’ 

x q y q - X q y q + X q -y q x q y 2 q - x q y q + x q - y q 
x q y q -y 2 + y q -l ’ x q ( -x q y q + y 2 -y q + l) 

Let us show how to find the analytic branch of log (x q —y q ) in the physical region. The values 
of the argument cover a region in the second, third and fourth quadrant of the complex 
plane, positive real axis included, intersecting the cut of the logarithm. An analytic branch 
can be found, for every polylogarithm, subtracting the discontinuity of the function across 
the branch cut. In the present case Disc(log(x)) = 2iri so that the analytic branch, log*(x), 
is 

log*(x q - y q ) = log (x q - y q ) - 2m 6(%(x q - y q ))), (D.3) 

where 0(x) is the Heaviside 9 function. Note that employing 6 functions we introduce a 
(removable) discontinuity on the real axis. In the next section we show how to remove this 
discontinuity. 

For Li 3 (aj) the cq’s lie on the branch cut whenever 9?(aj) > 1. We remove the cut 
ambiguity recalling that the physical region is defined via the Feynman prescription, and 
in order to get an analytic function we use the analytic continuation of Li 3 whenever 
9?(cii) > 1. This can be achieved, again, by means of 9 functions, 

Li 3 (ctj) = 6(1 - 3ft(a i ))Li 3 (a i ) + 6>($ft(a;) - l)C 3 (a,, a*), 
where cq is the sign of the imaginary part of a* due to Feynman prescription, 

<? = {-!,-1,-1,1,1}, (D.5) 


(D.4) 


(D.2) 
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and C 3 is the analytic continuation of IU 3 


C z{ai,(T.j) = Li 3 ( — ) - - log 3 (a*) + cq-m log 2 (a*) + )!-7r 2 log(aj). (D.6) 

\ai) 6 2 3 

D.l Analytic branches in the physical region 

The representations (D.3) and (D.4) are analytic everywhere in the physical region, but 
if the argument of the 9 function is zero the expression is not defined. Nevertheless, by 
construction, the new basis has no branch cuts and the limits to the discontinuous points 
are well defined. In order to remove these (spurious) discontinuities we impose the limiting 
values as long as the arguments of the 9 functions vanish. 

The limiting values to the real axis of Eq. (D.3) are such that log*(—1) = —in and 
log*(x > 0) is positive. This uniquely defines Eq. (D.3) over the entire physical region. 
With a similar reasoning we can find the analytic branch, log'(x), of the second logarithm 
of Eq. (D.l). Defining a = x q — y q and b = —y 5 x 2 + y q x q — y q + x q we arrive at 

log*(a) = [1 - 6 (9(a))] c(a) + 5 (9(a)) [3?(log(a)) - in], 

log f (&) = [1-S ($(&))] d(b) + 6 (3?(&)) [3?(log(6)) - in(l + 9 l/2 {x q , y q ))\ , (D.7) 


where 


and 


<5(x) 


1 forx = 0 
0 otherwise ’ 


(D. 8 ) 


c(a) = log (a) — 2in 6($s(a)), 

d(b) = log( 6 ) - 2 i 7 r 0(9(b)) - 2in6 1/2 (x q , y q )6(-$t(b)), 

Q 1 / 2 (x q , y q ) = 9o(n/2 - aig(x q ))9 0 (n/2 - arg (y q )), (D.9) 


where 0 < arg(x) < n. 

In Eq. (D.4), Li 3 (l) = C 3 (l,CTj), then 


Li^Oi) = 6 *o(l - »(o i ))Li 3 (a i ) + 0i(K(oi) - 1)C 3 (a i ,a i ), (D.10) 

where 0o(3 ; ) and 6\{x) are modified Heaviside 9 functions, such that 6 * 0 (0) = 0 and #i(0) = 1. 

Note that Eqs. (D.7) do not hold if x q , y q are both real, where the standard log has to 
be used. To avoid clutter we omit here the extra terms extending the validity of Eq. (D.7) 
to that region, while they are present in the ancillary files of the arXiv submission. Al¬ 
though the above representations might look cumbersome, they are nothing but simple 
combinations of elementary functions like the Heaviside 9 , and implementing them for nu¬ 
meric evaluations is straightforward computer algebra. The expressions in the ancillary 
files make only use of built-in Mathematica functions. 
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